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Abstract. A fermion realization of the compact symplectic sp{4) algebra provides 
a natural framework for studying isovector pairing correlations in nuclei. While these 
correlations manifest themselves most clearly in the binding energies of 0^ ground 
states, they also have a large effect on the energies of excited states, including 
especially excited 0+ states. In this article we consider non-deformed as well as 
deformed algebraic descriptions of pairing through the reductions of sp(q) (4) to different 
realizations of (2) for single-j and multi-j orbitals. The model yields a classification 
scheme for completely paired 0"'' states of even-even and odd-odd nuclei in the 
l'^3/25 1/7/21 ^nd I/5/22P1/22P3/2I59/2 shells. Phenomcnological non-deformed and 
deformed isospin-brcaking Hamiltonians arc expressed in terms of the generators of 
the dynamical symmetry groups 5*^(4) and Spq{A). These Hamiltonians are related to 
the most general microscopic pairing problem, including isovector pairing and isoscalar 
proton-neutron interaction along with non-linear interaction in the deformed extension. 
In both the non-deformed and deformed cases the eigenvalues of the Hamiltonian are fit 
to the relevant Coulomb corrected experimental 0^ energies and this, in turn, allows 
us to estimate the interaction strength parameters, to investigate isovector-pairing 
properties and symmetries breaking, and to predict the corresponding energies. While 
the non-deformed theory yields results that are comparable to other theories for light 
nuclei, the deformed extension, which takes into account higher-order interactions 
between the particles, gives a better fit to the data. The multi-shell applications 
of the model provide for reasonable predictions of energies of exotic nuclei. 
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1. Introduction 

The pairing problem, which was introduced first in atomic physics P, was later applied 
to nuclear physics [21 E] in an attempt to describe binding energies of nuclei and their low- 
lying vibrational spectra IHIHl- Recently, there has been renewed interest in this problem 
because of new experimental studies of exotic nuclei with relatively large proton excess 
or with N Z. This revival of interest in pairing follows from the recent development of 
radioactive beam facilities and attempts to bridge from nuclear structure considerations 
to astrophysical phenomena [HI E] • 

Along with approximate mean field solutions (for a review see jH]), the pairing 
problem can be solved exactly by means of various group theoretical methods, which 
allow one to explore the underlying symmetries. The SU (2) seniority model CHI 
provides for a good description of nuclei with large proton or neutron excess, where the 
like-particle pairing plays a dominant role. The simple "quasi-spin" {SU{2) ~ Sp{2)) 
approach ^2] not only offers an elegant way to understand the results from the 
conventional seniority scheme [H [121, based on f/(2j + 1) D Sp{2j + 1) D 5*0(3), 
but allows for a straightforward expansion to 5*0(5) ~ 5*^(4) to include protons and 
neutrons, which otherwise has proven to be too complicated. The generalization to the 
5*0(5) model O [T^l [IHl [HI introduces a relation between identical-particle and proton- 
neutron (pn) isovector (isospin r = 1) pairing modes. The addition of an isoscalar 
(r = 0) pn pairing channel is described within the framework of the 5*0(8) model 
[IHl[ini[2I]] and the Interacting Boson Model (IBM) [25. 

In the limit of dominant isovector pn pairing correlations, a simple 5*0(5) seniority 
model [221 122] is suitable. Our goal is to investigate properties of the isovector pairing 
interaction within the context of a fermion realization of the symplectic sp{4) algebra 
[which is isomorphic to so(5)]. The model space consists of J'^ = states with pairs 
coupled to isospin r = 1; mixing with t = pn pairs is not included. The importance 
of the isovector pairing for binding energies is suggested by experimental data, namely 
a r = 1 ground state for most N = Z odd-odd nuclei with mass number A > 40 
[211 [23 [23, and by the results of various theoretical studies [2Zll2Hll2ni[Sni[nil[221[221[211- 

While coupling to the isoscalar pn pairing mode may be important in some cases 
[23 ISni, we exclude - to the best of our ability - the ground states of nuclei that show 
fingerprints of isoscalar pairing correlations. In short, the Sp{4) model is applied to 0''" 
ground states of even-A nuclei and to the higher-lying O"*" isobaric analog states in most 
of odd-odd nuclei. We refer to these states as isovector-paired states. In this regard, it is 
important to note that the two-body interaction includes an isoscalar term in addition 
to the dominant isovector pairing interaction. The isoscalar pn force is related to the 
symmetry energy and is diagonal in the isovector-paired basis states with good isospin 
[2Z1. Diagonal high- J components of the nuclear interaction are also present in the 
model. The shell structure and its dimension play an important role in the construction 
of the fermion pairs and their interaction in accordance with the Pauli principle. The 
isovector pairing term is assumed to be particle-hole symmetric, which enters naturally 
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in points to a decrease in energy with respect to the mean-field solution of a no-pair 
theory [SH] . 

Limiting cases of sp{4) correspond to different reductions of sp{A) to u(2) and 
show distinct properties of different coupling modes of the isovector pairing interaction; 
specifically, proton-neutron [pn) and like-particle [pp and nn) pairing phenomena. The 
theory provides for a classification of nuclear states with respect to the number valence 
protons and neutrons occupying a major shell. The notion of a dynamical symmetry 
extends this picture to include an isospin-breaking phenomenological interaction which 
is related to a general microscopic Hamiltonian for the pairing problem. The final result 
can be written in terms of the second order invariants of the subalgebras of sp{4) which 
then reduces the problem to an exactly solvable theory. 

A g-deformation of the classical algebraic structure is introduced. A quantum 
extension of the dynamical symmetry approach is realized leading also to an exact q- 
deformed solution of the problem and its limiting cases. The motivation behind the 
g-deformed generalization is that, as a novel and richer model, it allows us to include 
non-linear features of the interaction and to investigate the respective changes this may 
require in the strength parameters and in the pairing gaps. Existing applications of the 
g-deformed algebraic structures to the pairing problem j^lj are restricted mainly to the 
SUq{2) limit [101 HI] of the dynamical symmetry approach presented here for Spqi^A). 

Fully-paired even-even and odd-odd nuclei, 32 < A < 100, in 1^3/2, I/7/2 and 
1/5/22^1/2 2^3/21^79/2 orbitals are considered in details in this investigation. An analysis 
of the results, obtained by fitting model parameters to experimental data both in the 
deformed and non-deformed cases, provides for a reasonable prediction of the relevant 
O"*" state energies of nuclei classified as belonging to a major shell and gives insight 
into their pair structure and isospin mixing. It also estimates the broad limits of 
applicability of such a simple algebraic model and its deformed non-linear extension- 
which is in agreement with the results of other theoretical approaches for describing 
pairing phenomena in nuclear systems [131 122 EH 1121 1121 Hlj • 

The paper is organized as follows. In the next section, the algebraic structure of the 
fermion realization of sp(4) and its deformation spg(4) is introduced with emphasis on 
the physical interpretation of the generators. In Section 3 the application of the algebraic 
constructions is realized through the introduction of a model Hamiltonian in both the 
deformed and non-deformed cases. In Section 4, the parameters of the Hamiltonians are 
presented as output of a fitting procedure to the respective experimental energies and 
the results are analyzed. A summary of our findings and the main conclusions as well 
as possible further developments of the approach are discussed in the final section. 

2. Algebraic structure of non-deformed and q-deformed sp(4) 

To introduce notation, we start with a brief review of the algebraic structures that 
enter into the discussion jU]. The sp(4) algebra is realized in terms of the creation 
(annihilation) fermion operator c]^^^ (cm,,^) cr = ±1, where these operators create 
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(annihilate) a particle of type cr in a state of total angular momentum j = i, |, |, 
with projection m along the z axis {—j < m < j). They satisfy Fermi anticommutation 
relations 

{Cm',<7', 4^,^} = Sm',raSa',cT A ^ 4n,a} = {Cm' , Cm,a} = (l) 

and Hermitian conjugation is given by (cj^o-)* = Cm,a- For a given a, the dimension of 
the fermion space is 2Qj = 2j + 1. 

The deformation of the spq{4) algebra is introduced in terms of g-deformed creation 
and annihilation operators aj^^. and a^.o-, (o^m.o-)* ~ (^m,a, where 
the limit g — ^ 1. The deformed single-particle operators are defined through their 
anticommutation relation for every a and m [IH]: 

t ±— t 

{«m,a> «m',a'} = {"m.a, 0!m' = 0, 

where by definition the g-anticommutator is given as {A, B}^k = AB + q^BA. A prop- 
erty with physics impact is the dependence of the deformed anticommutation rela- 
tions on the shell dimension and the operators that count the number of particles, 

^±1 = E^m=-iCL,±lCm,±l. 

Generators of the symplectic group — The generators of 5*^(4) and Spq{'i) are 
expressed in terms of non-deformed ^3 ESI ^^"^ deformed single-particle operators 
respectively, 
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where a, a' = ±1 and v4o,±i = (-Bo,±l)^ -^o,±i = (Co.ii)^- These operators create (anni- 
hilate) a pair of fermions coupled to total angular momentum and parity = 0"^ [T^ E] 
and thus constitute boson-like objects. The rest of the generators of 5*^(4) are D^ ^,.' = 

and for 5*^5(4) they are -E±i,^i = -7= I]l=-j "ln,±i"m,q=i, in 



addition to the number operators, N±i, which remain non-deformed in this realization 
of spq(4). The ten non-deformed (deformed) generators close on the symplectic sp(q)(4) 
algebra with the commutation relations given in |3SI- For nuclear structure applications 
we use the set of the g-deformed commutation relations that is symmetric with respect 
to the exchange of the deformation parameter q <-> q~^. 



Physical interpretation of the generators — When considered to be a 
dynamical symmetry, the 5^(4) symplectic group can be used to describe distinct 
collective nuclear phenomena through different interpretations of the a quantum number 
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in the 5*^(4) D ® SU{2) reduction. When a is used to distinguish between 

protons (cr = 1) and neutrons (cr = —1), the Cartan generators of the 5'p(4) group 
A'^-ti (with eigenvalues N±) enter as the number of the valence protons and valence 
neutrons, respectively. 

The significant reduction limits of sp(g)(4) are summarized for the non-deformed 
and g-deformed cases in Table 1, where by definition [X]k = gfclg-fc ^^'^ P± = 
(q"^^ + g^2^)/2. Table 1 consists of four different realizations of a two-dimensional 
unitary subalgebra w^q)(2) D u^{l) © sm^^^(2) (/i = {r, 0, ±}) and the corresponding 
second-order Casimir invariant of suq{2). The "classical" formulae are restored in the 
limit when q goes to 1. In the first realization, su^(2), the generators ro,± are associated 
with the components of the isospin of the valence particles. The SU^{2) limit describes 
proton and neutron pairs (pn), while the SU^(2) limit is related to coupling between 
identical particles, proton-proton (pp) and neutron-neutron (nn) pairs. 



Table 1. Realizations of the unitary subalgebras of sp(g)(4): isospin symmetry [fi = t), 
pn [ii = 0) and pp{nn) (/i = ±) coupling, along with the Casimir invariants of su^{2). 
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Within a representation fij, the space of fully-paired states is constructed by the 
pair-creation g-deformed operators -Fo,±i (jll) (non-deformed operators Ao,±i (13), acting 
on the vacuum state [IH] : 

Im, no, n_i)^ = {F^r^ {For (^^-i)""^ |0) , (5) 



where ni,?7.o,n_i are the total number of pairs of each kind, pp, pn, nn, respectively. 
The basis is obtained by orthonormalization of (0). The g-deformed states are in general 
different from the classical ones and coincide with them in the limit g — 1. 

The generalization of the pairing problem to multi-shells dimension El El leads 
to a natural expansion of the fermion realization of the sp(4) algebra, allowing the 
nucleons to occupy a space of several orbits. The commutation relations between the 
ten non-deformed (deformed) generators of the generalized Sp(^q^{4) and the related 
algebraic formulae (derived in the single- level realization [IH]) remain the same with the 
substitution Qj Q, where 2Q = 2Qj = '^j{2j + 1). 

3. Theoretical model with 5*^(4) dynamical symmetry 

In the deformed and non-deformed cases, the basis states \ni,no,n^i)^^^^ (0) give the 
isovector-paired 0"^ states of a nucleus with = 2ni + no valence protons and 
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= 2n„i + no valence neutrons. This yields a simultaneous classification of the nuclei in a 
given major shell and of their corresponding isovector-paired states. The classification 
scheme is illustrated for the simple cases of 1^3/2 with ^lj=3/2 = 2 (Table 2a) and I/7/2 
with flj=Y/2 = 4 (Table 2b). The total number of the valence particles, n = + A^„, 
enumerates the rows and the eigenvalue i of the third projection of the valence isospin tq 
enumerates the columns. Isotopes of an element are situated along the right diagonals, 
isotones - along the left diagonals, and the rows consist of isobars for a given mass 
number. The shape of the table is symmetric with respect to i (with the exchange 
rii n_i), as well as with respect to n — 2Q (middle of the shell). This is a consequence 
of the charge independent nature of the interaction and the Pauli principle, respectively. 

Table2a. Classification scheme of nuclei, ^3/2 = 2. 
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Table2b. Classification scheme of nuclei, ^7/2 = 4. The shape of the table is symmetric 
with respect to the sign of i and n — 2Q. The basis states for each nucleus are labeled 
by the numbers of particle pairs \ni,no,n^i). 
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Model Hamiltonian — As a natural approach within a microscopic picture, the 
most general Hamiltonian of a system with 5*^(4) symmetry, which preserves the total 
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number of particles, can be expressed through the group generators as following |46j : 

- - Z,(.„.o - (6) 

where G, F, E, C and D are phenomenological constant interaction strength parameters 
(G > 0, -F > for attraction), e > is a Fermi level energy. 

An important feature of the phenomenological Hamiltonian (jUj) is that it not only 
breaks the isospin symmetry {D ^) but it also mixes states with definite isospin 
values {F G). This is different from other applications of non-deformed and deformed 
sp(4) or o(5) algebras with isospin-invariant Hamiltonians ^1I221I1Z1- Although the 
degree of mixing is expected to be smaller than for isoscalar-isovector mixing, it may 
still add an interesting contribution to the study of the isospin mixing 13211^01 • 

Possible applications of the Hamiltonian to real nuclei can be determined through a 
detailed investigation of the various terms introduced in (jH)). The first two terms {G, F) 
of the Hamiltonian © account for J = isovector pairing between non-identical and 
identical particles, respectively. To refiect the assumption that a zero pairing energy 
corresponds to a state with no possible breaking of a pair jHH], a particle-hole concept is 
incorporated in these two terms (but not in the e-, G- and D-terms). Hole pair-creation 
(annihilation) operators can be introduced not only for identical particle pairs {pp or 
nn) j2Hl; but also for pn pairs. This corresponds to a change from the particle to the 
hole number operator, N± 2Q - N± for N± > Q and N ^ 4Q - N for N > 2Q. 

The next term (E) can be related to the symmetry energy P^IT^ as its expectation 
value in states with definite isospin is 

I f K..>) ^ (e-^) ^ E^JZ^, (7) 

which enters as a symmetry term in many nuclear mass relationships |^ 1^ . The 
second order Casimir invariant of sp(4) jSSl sets linear dependence between the terms in 
(jni), which yields to a direct relation between the symmetry and pairing contributions: 
a fact that has been already pointed out in a phenomenological analysis based on the 
experimental nuclear masses and excitation energies [S3]- 

The two-body interaction in ©, which is written in terms of the group generators, 
arises naturally from the microscopic picture. In the single-j case its form is jHj 



2 

{a} M,m,m' 



m',(T3 '--jr,m',o'4 1, 



where {a} = {(cti, as, 0-3, 0-4)} = {(+,+,+,+), (+,-,+,-), (+,-,-,+), 
(—,—,—,—)}. The coefficient (72 | V| (J4, (T3) is the expectation value of the two- 
body interaction potential between pairs of quantum numbers (T4, (T3 and ai, o"2. The 
second sum in ((HJ can be expanded into three terms. The first term corresponds 
to pairing to total angular momentum J = (M = 0), the second term includes 
high-J (J 7^ 0, m = m') components of the interaction and can be represented 
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by {r+,r„} - N/{2n), N{N - 1) and ToTq - N/i. The rest of the sum is the 
residual interaction that is neglected. A multi-shell generalization of the microscopic 
Hamiltonian (jH)) can be related to the phenomenological one © and the interaction 
strengths can be obtained in terms of the phenomenological parameters 

J = pairing / (+ + l^^'ol + +) = (-- I^PoI - -) = F/n 
interaction 1 (- + |VpJ + -) = (- + |VpJ -+) = G/fi 

(+ + \V\ + +) = {--\V\--) = C + D/2 

(- + \V\ -+) = 2C-D, (- + \V\ + -)= E/Q. 

This connection with the interaction matrix elements gives a real physical meaning 
to the constant phenomenological strength parameters, and, therefore, their estimation 
can lead to a microscopic description of the nuclear interaction. 

The three terms C, D and E in ^ that arise from the dynamical >S'p(4) symmetry 
are related to the microscopic nature of the pn isoscalar correlations. As can be clearly 
seen from the expression A(_^2 + M + lM^)_(c+ A)M^_P_^)(^2_^) [g^^ 

(jH))] and from relation Q, for D = E /2VL and C + D /2 = Owe obtain the J-independent 
pn isoscalar force. It is closely related to the symmetry energy (E), is diagonal in the 
pairing basis with G = F and can be compared to j^ZlEll- Therefore, the S'p(4) model 
interaction consists of isovector [pp, nn, pn) pairing and isoscalar (pn) force in addition 
to a possible isospin-breaking term and J > identical-particle pairing correlations. 

In this way, the phenomenological Hamiltonian © can be used to describe general 
properties of the nuclear interaction, which serves as a motivation to fit the theoretical 
expectation values of (jHl) to the energies of the corresponding 0"^ states of nuclei in a 
very broad region. 

Within the algebraic framework, the important reduction chains of the symplectic 
algebra to the unitary two-dimensional subalgebras allow the Hamiltonian (jH)) to be 
expressed through second-order operators CJ'*^'^ (Table 1): 

H = - ri^Cl - ri^rl - r,^Cl - r^.D^ - v^iC^+C^) - V6{D+,{D+,-l) + 
D^iiD^i-l))-r]jN + r]8. (10) 

The ?7j-coefficients (i = 1, 2, . . . , 8) in ()10|) are not linearly independent; they are related 
to the phenomenological parameters of the model (jHl) in the following way: 



V7 



■ C(l - 4fi)/2 - D/A -{E- G)/(4fi), N <2Vt 
C{1-AVl)/2-D/A+{E-G)/{AVl), N>2VL] 



2Cn'^ + G/2 , N<2n 

2Gn^ + G/2 + {E -G), N >2n. ^ ' 

The ratios ?72/'7i, V-i/vs^ Ve/V5 determine the extent to which the symmetry in each limit 
is broken 
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In the g-deformed Hamiltonian can be constructed that is analogous to (|l(Jj) 

and is chosen to coincide with the non-deformed one © in the limit q ^ 1 



E,{{T+,T^} 



-Cq2fl 



2f2 



' 20. 



N 
(12) 



where e'' = e'' + — 2Vt)Cq + ^ > is the Fermi level of the nuclear system, Kq is 
related to (Table 1), Gg, Fg, Eg, Cq and Dq are constant interaction strength pa- 
rameters and in general they may be different than the non-deformed phenomenological 
parameters. 



Matrix elements of the Hamiltonian — In the SU^{2) limit (pn-coupling) the 
energy eigenvalue of the non-deformed pairing interaction GAqBq is 
G 2Q-n + no + l G 



[n 



2z/o)(4fi-ra-2z/o + 2) 



(13) 



2 8fi 

and in the SU'^{2) limit (like-particle coupling) the energy of the non-deformed pairing 
interaction FA±iB±i is 

F F 
Sppinn) = ^^±i(^ + n±i-N± + l) = — {N± - i^i) (2fi -N±-ui + 2). (14) 

In each limit, uq = ni + n^i and ui = no are the respective seniority quantum numbers 
that count the number of remaining pairs that can be formed after coupling the fermions 
in the primary pairing mode and they vary by Ai^o,i = 2. 

To investigate the influence of the deformation on the pairing interaction, the 
eigenvalue of the deformed pairing Hamiltonian is expanded in orders of x (g = e'^) 
in each limit 
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(16) 



where the non-deformed energies (fT!?|) and (fT^ are the zeroth order approximation of the 
corresponding deformed pairing energies. While the proton-neutron interaction is even 
with respect to the deformation parameter x, the identical particle pairing includes odd 
terms as well through the coefficient p±. The expansions in the pairing limits fpSj) and 
(|TH|l ) introduce non-linear terms with respect to the pair numbers, space dimension and 
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the non-deformed pairing energies. They serve as a simple example of the contribution 
of the g-deformation compared to the non-deformed model, which is a straightforward 
result of the quantum definition. 

In general, the Hamiltonian (jH)) is not diagonal in the basis set (Table 2b). The 
linear combinations of the basis states describe the spectrum of the isovector-paired O"*" 
states for a given nucleus. The pairing Hamiltonian Hpair (© with E = C = D = 
and e = 0) gives a transition between the states with different kinds of pairing while 
preserving the total number of pairs, A^, that is, two pn pairs scatter into a pp and a nn 
pair, and vice versa 

G 

\Hpair\ \ni,no,n_i) = {Spn + Epp + Enn) \ni,no,n_i) - -^ni^^i 1^1 - l,no + 2,n_i - 1) 

F 

- Q^oi^o - 1) |ni + 1, no - 2, n_i + 1) , (17) 

where Epn^pp^nn are given in (fT^ and (fT^ and ni,no,n_i are particle or hole pairs. 
We are also able to find an analytical form of the g-deformed analog of (fTTjl 

G 

\Hg^pair\ 1^1, no, n_i) = (5^„ + E^p + eU |ni, no, n_i) - -^"-i^-i I'^i - 1, rio + 2, n_i - 1) 

P / "0 — 1 

- Q^^f^ *^''(^) + I'^o - 2,n_i + 1) , (18) 
^ -I fc=i 

where E^^ pp^^^ are given in and jTH)). We define n±i = pj^p+p_ t^^^^] [^"±1-^-1/2] _i_ 
n±„ S,{k) ^ [2,_n-i/2]^ Eto ¥[2,_i_,]^ 4A: and [2^]^ ^ ^ 2. 

211 

4. Applications to nuclear structures 

Jf-.l. 0^ -state energy for even-A nuclei. Discussion of the results. 

The eigenvalues of the Hamiltonians © and (fT^ describe nuclear isovector-paired O"*" 
state energies, which are fit to experimental values jSHllSZl- For even-even nuclei and 
for some odd-odd nuclei {Z ^ A^), the lowest O"*" state is the nuclear ground state and 
the positive value of its energy is defined as the binding energy, \BE\. The binding 
energy of a nucleus is an important quantity because it is related to the nuclear mass 
and lifetime. Other odd-odd nuclei have a higher-lying O"*" excited state which is an 
isobaric analog of the corresponding even-even neighbors. 

The phenomenological parameters in © and (fT^ are determined by a non-linear 
least-squares fit of the lowest isovector-paired O"*" state energies (maximum eigenvalues 
of I if I © and \Hq\ (O) to the Coulomb corrected experimental values: 

Eo^^P(iV+, iV_) = |Eexp(iV+, N_)\ - l^explcore + Vcoul{N+, N_), (19) 

where the binding energy of the core l-Eexpl core is subtracted in order to focus only on the 
contribution from the valence shell. The energies need to be corrected for the Coulomb 
repulsion since it is not accounted for by the model Hamiltonian. In (fT^ the Coulomb 
potential is taken relative to the core and is derived in |HH] . 
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The parameters and statistics, obtained from the fitting procedure, are shown in 
Table 3. In both the non-deformed ("non-def column) and deformed cases ("g-def 
column), three groups of even- A nuclei are considered: (I) 1*^3/2 = 2) with a core f^S 
(Table 2a); (II) I/7/2 {^l = 4) with a core |°Ca (Table 2b); and (III) 1/5/22^1/22^3/2 lc/9/2 
major shell [Q = 11) with a core 28^^- Iii each group, the number of the valence protons 
(neutrons) varies in the range N± = 0, . . . , 2Q and the total number of nuclei that enter 
into the Sp^i) systematics is 2fi(fi + 1) + 1 (13 for (I), 41 for (II) and 265 for (III)). 
The residual sum of squares S = (|-E'o''| — |-E'o^''|)^ and the chi-statistics x = \J ^f-n 
define the goodness of the fit, where rip is the number of the fitting parameters and A'^^ 
is the number of nuclei with available data (A^^^ is 13 in (I), 36 in (II) and 100 in (III)). 

Analysis of the results (Table 3) shows that for 1^3/2 the pairing parameters are 
almost equal {G ~ F) as it is expected for light nuclei, and they differ, G > F, for 
1/7/2 by 0.07 and for (III) by 0.06. Based on the estimation of the parameters (Table 
3) and the correlations (|TT| the extent to which the symmetry in each limit is broken 
can be evaluated. In the limit SU'^{2), the breaking of the isospin invariance 772/771 
is in general small for light nuclei {ri2/r]i = 0.090 for 1^3/2 and ?72/''7i = 0. 133 for 
1/7/2); which is in agreement with the experimental data for this region. For medium 
nuclei in the l/5/22pi/22p3/2l5'9/2 major shell the isospin breaking is significantly greater, 
V2/V1 = 0. 628. 

Table 3. Fit parameters and statistics. G, F, C, D, e and x are in MeV, S is in MeV"^. 
Quantities marked with the symbol * are fixed for a given fit. 





(I) 


(11) 


(III) 




1(^3/2 


1/7/2 


1/5/22^1/22^3/2 16-9/2 




non-def 


g-def 


non-def 


g-def 


non-def 


g-def 


X 


0* 


-0.015 


0* 


0.124 


0* 


0.215 


g = 


1* 


0.985 


1* 


1.132 


1* 


1.240 


G/Q 


0.709 


0.709* 


0.525 


0.525* 


0.352 


0.352* 


F/Q 


0.702 


0.702* 


0.453 


0.453* 


0.296 


0.296* 


C 


0.815 


0.815* 


0.473 


0.473* 


0.190 


0.190* 


D 


-1.282 


-1.282* 


-0.971 


-0.971* 


-0.796 


-0.796* 


E/{2VL) 


-1.409 


-1.409* 


-1.120 


-1.120* 


-0.489 


-0.489* 


e 


9.012 


9.012* 


9.359 


9.359* 


9.567 


9.567* 


S 


1.720 


1.719 


16.095 


15.673 


300.284 


238.280 


X 


0.496 


0.378 


0.732 


0.669 


1.787 


1.551 



An observation about the pn-pairing strength is that most of pn-coupling study has 
been done assuming good isospin, that is F = G. However, the free nucleon-nucleon 
data jnni indicates that the r = 1 pn-pairing strength (G) is slightly bigger than the 
like-particle pairing strength (F), which is confirmed by the fits presented in Table 3 
and is investigated in other studies j2Zl ED]- There are many different values for the 
like-particle pairing strength used in literature. The most common value is taken to be 
proportional to 1/A jSl CHI EH IBS EBj and is consistent with the experimental pairing 
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Figure 1. Coulomb corrected 0+ state energy Eq vs. the isospin projection i 
for tlic isotopes of nuclei with Z = 20 to Z = 28 in the I/7/2 level, $^7/2 = 4. 
The experimental binding energies E^^^ (symbol "x") are distinguished from the 
experimental energies of the isobaric analog 0+ excited states E'^'^ (symbol "o"). 
Each line connects theoretically predicted energies of an isobar sequence. The nuclei 
for which experimental data is not available are represented only by their predicted 0+ 
state energy. 



gaps derived from the odd-even mass differences [02 EH] • The values of F, obtained by 
our theoretical symplectic model, fall within the limits of their estimation. In this way, 
they are expected to reproduce the low-lying vibrational spectra of spherical nuclei in 
the SU^{2) limit. When the results from all the three non-deformed fits are considered 
(Table 3), the identical-particle parameter F/VL is found to decrease with the mass 
number as 23. 9/yl. In a similar way, one can find the dependence of the pn pairing 
strength parameter G/Vt on the mass number to be 25. 7 /A. 

The estimate of the parameters (Table 3) reveals the properties of the nuclear 
interaction as interpreted by connection (jH)). The J = pairing interaction (Vpg) is 
always attractive, while the overall high-J component of identical-nucleon coupling 
(± ± \V\ ± ±) might be repulsive. The J > proton-neutron "direct" interaction 
(- + \V\ - +) is attractive, but not the "exchange" part of it (- + \V\ + -)(£'< 0). 

In all cases there is a good agreement with experiment (small x) 1 can be seen in 
Table 3, as well as in Figure ^ for region (II). Part of our results, namely for the binding 
energies (but not for the excited 0^ state energies), can be compared to other theories. 
A direct comparison of the chi statistics is impossible because of the different data sets 
and energy levels determined by the various theories. However, if we select only the 
data subsets that are equivalent for the nuclei in the 1^3/2 and/or I/7/2, our results are 
much closer to the experimental numbers than those for the Hartree-Fock-Bogoliubov 
(HFB) model [IHj and the semi-empirical model jH] and comparable with those of the 
jj-coupling shell model ^21 112] and the isovector and isoscalar pairing plus quadrupole 
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model [SZ|- In this way the simple 5*^(4) model is tested and proves its validity when 
applied to light nuclei in single-j level. In this region many symmetries are conserved 
allowing for a possible reduction of the number of fitting parameters. However, the 
free parameters in the fits presented in Table 3 reflect the symmetries observed in light 
nuclei and the non-negligible symmetry breaking in medium-mass nuclei. 

Properties of the pairing interaction — A model with the 5*^(4) dynamical 
symmetry permits an independent investigation of the different kinds of pairing 
interactions in the limiting cases of the non-deformed p3|) . (fT^ as well as the deformed 
versions (fTH|l . (fTT)|l of the theory. In the S'f/^(2) limit, the symplectic model reproduces 
the properties of the identical- nucleon pairing {Spp + (fT^ . for which the usual 
parabolic dependence of £pp(nn) on N± holds [101 IHl UHl IHHj- The dependence of the 
like-particle energy on the isospin projection i (Figure El^a)) reveals another property of 
the pairing mode, a Ai = 1 staggering of the identical-nucleons pairing energies of the 
odd-odd and even-even nuclei. In contrast with this, the pn limit (ep„) shows a smooth 




Figure 2. Identical (G = 0, F = 1.81 (a)) and non-identical (G = 2.10, F = (b)) 
particle pairing energies vs. i for isobars with A ~ AO to A ~ 56 in the I/7/2 level 

behavior (Figure El (b)). The SU^{2) limiting case yields a proton-neutron coupling 
that has its maximum when = N_ [i = 0), which is consistent with a— clustering 
theories EZj and the charge independence in the region of light nuclei when protons 
and neutrons fill the same shell [221 In both limits {SU^{2) and SU^{2) ), the pairing 
energy decreases when the difference between proton and neutron numbers increases. 

In most nuclei the different pairing interactions coexist and the contribution of each 
of the pairing modes, {Hpp + Hnn) = F {A^iB^i + A_iB_i) and {Hpn) = G {AqBq), to 
the total pairing energy (Hpair) = (Hpn) + {Hpp + Hnn) can be investigated. A AN = 2 
staggering exists for both pairing interactions (Figure Efa)) [22linHlEI|- For A+ = 
odd-odd nuclei the t = 1 pn pairs give the dominant contribution, while for the even- 
even A_|_ = N_ nuclei both pairing modes contribute almost equally with a slightly 
greater like-particle contribution. Although difference between even-even and odd-odd 



An Algebraic Pairing Model with 5*^(4) Symmetry and its Deformation 14 

nuclei exists in each pairing contribution, the total pairing energy has a surprisingly 
smooth behavior. The contribution from the symmetry term [E term in ©) [23 El] 
restores the staggering, as it decreases the energy of the odd-odd nuclei with respect to 
their even-even neighbors (Figure El^b)). 

Rough measures for the number of pn and like-particle pairs are the quantities 
^ (Hpn) and {Hpp + Hnn), respectively, which are related to the pairing gaps [221 122] ■ 
The "number" of pn pairs (Figure |3)) is bigger than the "number" of pp{nn) pairs for 
odd-odd = nuclei, and is of the same order as for the even-even nuclei [22linHllSIl- 

Deformed non-linear model — To investigate the role of the g-deformation, 
we performed again the fitting procedures for the same regions ((I), (II) and (III)) but 
using the deformed Hamiltonian ()12|1 . For each group of nuclei, the outcome of a fit with 
all possible parameters {Gq, Fg, Eg, Cg, Dg, e'^ and q) indicates that the introduction 
of the g-deformation does not vary the rest of the parameters. Based on this result, we 
considered the deformation to be independent of the other parameters and varied only 
q in the fit (the rest of the parameters were kept fixed with values obtained from the 
non-deformed fit). The results are shown in the "g-def columns in Table 3. The fits 
with and without a deformation can be compared by using the residual sum of squares 
(S*), which is always smaller in the deformed case (Table 3). 

Although it stands in contrast with other g-deformed applications [101 the 
decoupling of the g-deformation from the interaction strengths is not an assumption but 
results from comparisons to experimental data over total of 149 nuclei. It implies that 
while leaving the strength of the two-body interactions unchanged, the g-deformation 
allows one to take into account, in a prescribed way, complicated dependence of the 
energy eigenvalues on the number of nucleons/pairs and space dimension that cannot 
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Figure 3. Lowest isovector-paired 0+ state energies vs. when = for tlie 
nuclei in I/7/2: (a) pairing energies: pn, + nn and total; (b) symmetry energy (E 
term in IBJ) and total pairing + symmetry energy 
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Figure 4. Non-identical (pn) and identical {pp or nn) pairing "numbers" vs. for 
the N+ =^ nuclei with Z = 20 to Z = 28 in I/7/2 

be reproduced by any two-body interaction (for example, see fll5|) and (fTBj) ). Moreover, 
similar terms are expected to arise from higher-order interactions between the particles. 
In this way the g-parameter introduces some non-linear residual interaction not present 
in the two-body Hamiltonian 0. 

The observed independence of the pairing strengths on the g-parameter, suggests 
that while the deformation does not change the strength to couple two particles, it can 
model many-pair effects and can influence the energy spectrum. As an illustration, 
in each of the dynamical limits we investigate the quantities i?p„, = e'^^/Spn and 
Rpp+nn = {^Ip + £nn) I i^pp + ^nn) that givc an additional contribution to the pairing 
energy in the deformed case (Figure Ej) (compare to the analytical expansion with 
respect to x of the energies, (|T5|l and |T6|) ). In the limit of pn-pairing, Rpn does not 
signiflcantly change when q is close to one and it decreases for all g 7^ 1. The ratio 
Rpp+nn increases (decreases) monotonically with q only for nuclei with a primary pp (nn) 
coupling. Even though both SU^{2) groups are complementary, the different behavior 
of the multiplication constants p± (Table 1) is responsible for different impact of the 
deformation in various isotopes. This accounts for the differences in the experimental 
data between mirror nuclei even after the Coulomb energy correction. In the limit of 
identical-particle coupling, when q increases from one (g > 1) neutron pairs are less 
bound and proton pairs give a larger pairing gap, and vice versa for g < 1. In this way, 
the deformation parameter can determine the degree to which the pp coupling differs 
from the nn coupling. 

The signiflcance of the higher-order terms that enter through the g-deformed 
theory can be estimated through a comparison with experiment. In general, the fltting 
procedures determine values for x (Table 3) that are small. The reason may be that 
while higher-order effects may be signiflcant in nuclei they probably cancel on average 
when the g-parameter is one and the same for all nuclei. However, in two of the cases, 
(II) and (III), it is of an order of magnitude greater than the estimation of other physical 
applications ([201 and references there) and for the shell I/7/2, our value (g = 1. 132) 
is comparable to the values obtained in a g-deformed like-particle seniority model |4Uj : 
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q = 1. 1585 for the neutron pairs and q = 1. 1924 for protons. For the nuclei in the 
multi-j shell our model yields a bigger g-parameter than for the lighter nuclei in single- j' 
shell (Table 3), where the small number of valence nucleons is not sufficient to build 
strong non-linear correlations. This suggests that the g-deformation is more significant 
for masses A > 56. 



0.95 








0.9 








0.85 


L ;'/V 














0.8 




— (1.1) 








(1,3), (3,1) 


""Sc, 


0.75 










; ■'/ 


(2,2) 


44-pi 


0.7 




- (2,4), (4,2) 






; / 


(3,3) 


46^ 




(0,4) 




(1 ,3) 




(2,2) 


44-pi 


(2,4) 


46-pi 


(3,1) 




(3,3) 


46^ 


(4,0) 





1 



1 



Figure 5. Ratios Rpn and Rpp+nn vs. q for several nuclei with a typical behavior in 
the 1/7/2 level 



4.2. Predicted energies 

The fitting procedure not only estimates the magnitude of the pairing strength and 
describes the type of the dominant coupling mode, it also can be used to predict nuclear 
energies that have not been measured. From the fit for the I/7/2 case the binding energy 
of the proton-rich "^^Ni nucleus is estimated to be 348.19 MeV, which is by 0.07% greater 
than the sophisticated semi-empirical estimate of Likewise, for the odd-odd nuclei 
that do not have measured energy spectra the theory can predict the energy of their 
lowest 0+ isobaric analog state: 358.75 MeV (^V), 359.49 MeV {^'^Mn), 357.56 MeV 
{^^Co), 394.16 MeV (^°Co). The predicted energies are calculated for q = 1.132 (Table 
3 (II)) as the fit with deformation has a smaller uncertainty compared to the non- 
deformed one. The Sp^q) model predicts the relevant 0+ state energies for additional 
165 even- A nuclei in the medium mass region (III). The binding energies for 25 of them 
are also calculated in |13]. For these even-even nuclei, we predict binding energies that 
on average are by 0.05% (non-deformed case) and by 0.008% (for q = 1.240) less than 
the semi-empirical approximation |44j . 
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5. Conclusion 

We constructed a model with a symplectic dynamical symmetry group 5'p(g)(4) D 
f/(l) eg) S'[/(q)(2) in the non-deformed limit as well as in the g-deformed generalization. 
A phenomenological Hamiltonian was written in terms of the generators of the group 
and this in turn was used to describe pairing correlations in nuclei. The relation of 
this approach to a general microscopic pairing Hamiltonian was obtained. The theory 
was tested by fitting calculated energies to the relevant experimental 0^ state energies 
for single-j levels, namely 1^3/2 and I/7/2, and for a multi-j 1/5/22^1/2 2^3/2 lfi'9/2 shell. 
In general, the fitting procedure yielded results that were in good agreement with the 
experiment. The theory predicted the lowest 0"*" isovector-paired state energy of nuclei 
with a deviation of at most 0.5% in the energy range considered. It was used to predict 
the binding energy for even-even nuclei and the lowest isovector-paired 0"*" state energy of 
odd-odd nuclei in the proton-rich region. The phenomenological pairing parameters and 
the strength of the pairing interaction were determined. In agreement with experimental 
analysis, breaking of the isospin invariance, 'r]2/vi 7^ 0; isospin mixing, G > F, is 
observed, except for light nuclei in the 1^3/2 level. 

The theoretical model with S'p(4) dynamical symmetry and its g-deformed version 
was used to investigate in greater detail the properties of the isovector pairing 
interaction. The study reveals that the pairing energy along with the symmetry energy 
are responsible for the experimentally observed staggering between even-even and odd- 
odd nuclear isovector-paired 0"'' state energies. Overall, the results show that the 
symplectic model can be used to provide a reasonable description of isovector-paired 
0^ states in nuclei, confirming in its limit results previously published for like-particle 
pairing correlations. When only the nuclei with 0^ ground states are considered, the non- 
deformed 5*^(4) model is comparable, for the region of light nuclei, to earlier theories. At 
the same time it gives some insight into the study of symmetry breaking and isovector 
pairing correlations, and it is based on a simple approach that is applicable in a broad 
region of the nuclear chart, including odd-odd and exotic nuclei. 

The g-deformed case gives the best overall results. The 5'pg(4) dynamical symmetry 
approach yielded a g-deformed exact solution and we were able to derive the g-deformed 
matrix elements of the interaction in a simple analytical form. In addition to the 
broken symmetries of the non-deformed model, the g-deformation breaks the symmetry 
between protons and neutrons, which again is small for light nuclei and consistent with 
experiment. The introduction of g leads to a decrease of the like-particle pairing gap for 
neutrons and an increase for protons as g increases from one. The g-parameter was found 
decoupled from the interaction strength parameters. This observation suggests that 
while the deformation does not influence the two-body interaction, it introduces higher- 
order interactions between the particles, which are neglected in most non-deformed 
models. The g-deformation is mass and shell dimension dependent and its effects are 
more significant in the medium mass region. In the present study, the g-parameter was 
found as high as 1.240 and is expected to be greater if its influence is not averaged over all 
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nuclei in a major shell. This suggests the need for a more elaborate investigation of the 
role of the g-deformation in each individual nucleus and the relation of the g-parameter 
to the underlying nuclear structure. 
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